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Chapter 1 



Dynamic localization in optical 
lattices 



The concept of dynamic localization goes back to an observation reported by Dunlap and 
Kenkre in 1986: The wave packet of a single particle moving on a single-band tight-binding 
lattice endowed with only nearest-neighbor couplings remains perpetually localized when 
driven by a spatially homogeneous ac force, provided the amplitude and the frequency of 
that force obey a certain condition pQ. When trying to overcome the limitations of the model, 
it is comparatively straightforward to deal with an arbitrary form of the dispersion relation 
- thus abandoning the nearest-neighbor approximation — and with arbitrary time-periodic 
forces, thus doing away with the restriction to purely sinusoidal driving [2J. But in any real 
lattice system an external time-periodic force will induce interband transitions, and it is by 
no means obvious whether dynamic localization can survive when these come into play. 

In this chapter we consider ultracold atoms in driven optical lattices, which provide par- 
ticularly attractive, experimentally well accessible examples of quantum particles in spatially 
periodic structures exposed to time-periodic forcing [U SJ [5] . Such systems are much cleaner, 
and more easy to control, than electrons in crystal lattices under the influence of ac electric 
fields, for which the original idea had been developed pp. With the help of results obtained 
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by numerical calculations we illustrate that such ultracold atoms in kHz-driven optical lat- 
tices exhibit dynamic localization in almost its purest form if the parameters are chosen 
judiciously, despite the potentially devastating presence of interband transitions. 

When viewing dynamic localization as resulting from a band collapse [61 [7] , far-reaching 
further possibilities emerge. Namely, the actual strength of deviations from exact spatial 
periodicity, be they isolated [8], random [9], or quasiperiodic [TUl ITT] , is measured relative 
to the effective band width. Thus, when the band in question almost collapses in response 
to time-periodic driving, the effects of even slight deviations from exact lattice periodicity 
are strongly enhanced. This allows one, in particular, to coherently control the "metal- 
insulator" -like incommensurability transition occurring in sufficiently deep quasiperiodic op- 
tical lattices [TUl EJ 12] . While the very transition has already been observed with Bose- 
Einstein condensates in bichromatic optical potentials [13] , its coherent control by means of 
time-periodic forcing still awaits its experimental verification. 

1.1 The basic idea 

The one-dimensional tight-binding system described by the Hamiltonian 



where \£) denotes a Wannier state localized at the £th lattice site, and J is the hopping 
matrix element connecting neighboring sites, is about the simplest model for the formation 
of Bloch bands. Assuming that the unspecified number of sites is so large that finite-size 
effects may be neglected, its energy eigenstates are Bloch waves 



i 

labeled by a wave number k; the lattice period is given by a. The corresponding energy 
dispersion relation reads 
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here we assume J > 0, so that its minimum is located at k = mod 2ir/a. Now we let an 
external time-dependent, spatially homogeneous force F(t) act on the system, such that the 
total Hamiltonian becomes 

H(t)=H + H 1 (t) (1.4) 

with 

H^t) = -F(t)J2\£)a£(£\ . (1.5) 

i 

It is easy to verify that the wave functions 

= ex p(-^ jl&T E(q k (r))} J2 \ £ ) ex P Q*Qk(t)a) (1.6) 

then are solutions to the time-dependent Schrodinger equation, provided the time-dependent 
wave numbers qk{t) introduced here obey the "semiclassical" relation 

hq k (t) = F(t) . (1.7) 

We demand that quit) be equal to k at time t = 0, and therefore set 

1 /■* 



g fc (t) = A: + - / drF(r) . (1. 



These wave functions fll.6p . originally considered by Houston in the context of crystal elec- 
trons exposed to a uniform electric field superimposed on a periodic lattice potential [H], 
are known as "accelerated Bloch waves" , or Houston states. 

In the particular case of a monochromatic force with angular frequency u and ampli- 
tude F\, given by 

F(t) = Fx cos(ut) , (1.9) 



one has 



q k (t) = k + ^sin(ujt) , (1.10) 
huj 



so that qk{t) naturally acquires the temporal period T = 2tt/u of the driving force. Then 
also E(qk(t)) is T-periodic, but the Houston state (11.61) is not, because the integral appearing 
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in the exponential prefactor acquires a contribution which grows linearly with time. In order 
to extract that contribution, we calculate the one-cycle average 

e(k) = ^ £ dr E{q k {r)) 

= — 2 J cS cos(ka) , (1-11) 

thus obtaining an effective hopping matrix element given by 

J cS = Jh(^) , f 1.121 



with J (z) denoting the Bessel function of zero order. We then write 

exp(-^J*dTE{q k (T))^ =exp(-ijT t dr [E{q k {r)) - e{k)fj exp(- ie{k)t/h) , (1.13) 

so that the first exponential on the right hand side now is T-periodic by construction. Hence, 
for the T-periodic force ( II .9p the Houston states (I1.6P can be cast into a form 



\Mt)) = \u k (t)) exp( - ie{k)t/h) (1.14) 
with T-periodic functions \u k (t)), 

\u k (t)) = \u k (t + T)) . (1.15) 

This leads to a remarkable conclusion. Any wave packet governed by the full Hamilto- 
nian (jl.4p with periodic forcing (11. 9p can be expanded with respect to these states (11.141) 



with coefficents that are constant in time, because the time-dependence already is fully in- 
corporated into the states themselves. After each cycle T the T-periodic functions \u k (t)) 
are restored, so that the time evolution of the wave packet, when viewed stroboscopically at 
intervals T, is determined by the different "speed" of rotation of the complex phase factors 
exp(— is(k)t/h) of the packet's individual components. But if all quantities e{k) are equal, 
which according to Eqs. (II. lip and (11.121) occurs when the scaled driving amplitude 

K ^ (1.1.) 
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equals a zero of the Bessel function Jo, all phase factors evolve at the same speed, so that the 
wave packet reproduces itself exactly after each period: There is some T-periodic wiggling, 
but no long-term motion. This, in short, is dynamic localization pp. 

The above argument appears so special, and the decisive step f ll . 13j) so swift, that it is not 
easy to see how to transfer this finding to more realistic situations: How can one incorporate 
deviations from exact lattice periodicity into this reasoning? How to proceed when several 
bands are coupled by interband transitions? The answer to these questions is provided by 
the Floquet picture, which does not directly take recourse to the spatial lattice periodicity, 
but rather builds on the temporal periodicity of the Hamiltonian: When H(t) = H(t + T), 
there exists a complete set of solutions to the time-dependent Schrodinger equation of the 
particular form 



where the functions \u n (t)) = \u n (t + T)) inherit the T-periodicity of the underlying Hamil- 
tonian. These states are known as Floquet states; the quantities e n are dubbed as quasiener- 
gies [HI [IH [T7J [18]. Obviously the Houston states fjl.Sp with time-periodic forcing are 
particular examples of such Floquet states; from now on we employ an abstract state label n 
instead of the wave number k in order to also admit settings without lattice periodicity. 
In the case of the Houston-Floquet states, the determination of their quasienergies (11.111) 
essentially was a by-product of the solution of an initial value problem. The general case, 
however, has to proceed along a more sophisticated route: Floquet states and quasienergies 
are determined by solving the eigenvalue problem 



posed in an extended Hilbert space of T-periodic functions; in that space time plays the role 
of a coordinate. Therefore, if (ui(t)\u2(t)) is the scalar product of two T-periodic functions 
in the usual physical Hilbert space, their scalar product in the extended space reads [18] 





(1.17) 








Hence, we write \u n (t)) for a Floquet eigenfunction when viewed in the physical Hilbert 
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space, and \u n (t))) when that same function is regarded as an element of the extended space. 

A most important consequence of this formalism stems from the fact that when \u n (t))) 
is a solution to the problem fl 1 . 1 8 [) with eigenvalue e n , then \u n (t) exp(imut))) is a further 
solution with eigenvalue e n + mhu, where we have set u = 27r/T, and m is any (positive, zero, 
or negative) integer, in order to comply with the required T-periodic boundary condition. 
For m/fl these two solutions are orthogonal with respect to the scalar product fll . 19j) . But 
when going back to the physical Hilbert space, one has 

\u n (t) exp(imujt)} exp( — i(e n + m%u)/fij = \u n (t)) exp(—is n t/h) , (1.20) 

so that the two different solutions represent the same Floquet state fll. 171) . We conclude 
that a physical Floquet state does not simply correspond to an individual solution to the 
eigenvalue problem (|1.18p . but rather to a whole class of such solutions labeled by the state 
index n, whereas the "photon" index m distinguishes different representatives of such a class. 
Likewise, a quasienergy should not be regarded as a single eigenvalue, but rather as a set 
{e n + mhbj | m = 0,±1,±2, ...} associated with one particular state n, while m ranges 
through all integers. Therefore, each "quasienergy Brillouin zone" of width hu contains one 
quasienergy representative of each state. 

The time evolution of any wave function can then be written as a Floquet-state expansion, 

= E Cn\u n (t)) exp(-ie n t/h) , (1.21) 

n 

where the coefficients c n remain constant in time. This is one of the main benefits offered 
by the Floquet picture, and allows one to draw many parallels to the evolution of systems 
governed by a time-independent Hamiltonian. 

Equipped with this set of tools, it is now clear how to investigate the possible occur- 
rence of dynamic localization in realistic lattice structures: One has to solve the eigenvalue 
problem (11.181) for the Hamiltonian with the respective full lattice potential, and to enquire 
whether the resulting quasienergy bands collapse at least approximately, that is, acquire 
negligible widths for certain parameters. If so, any wave packet prepared in a quasienergy 
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band at a collapse point will suffer from "prohibited dephasing" , as in the archetypal model 
specified by Eqs. (II. ip . (II. 5p . and (jl.9P ; and thus remain dynamically localized. Interband 
transitions then are automatically included, with multiphoton-like resonances manifesting 
themselves through quasienergy-curve anticrossings [19] . 

In the following section we will carry through this program for ultracold atoms in driven 
one- dimensional optical lattices. 



1.2 Does it work? 



A one-dimensional optical lattice is created by two counterpropagating laser beams with 
wave number k^, suitably detuned from a dipole-allowed transition of the atomic species 
moving in this standing light wave. By means of the ac Stark effect, the spatially periodic 
electric field experienced by the atoms then translates into a cosine potential 

Vi a t(x) = yCos(2A; L x) (1.22) 

for their translational motion along the lattice, with a depth Vo that is proportional to the 
laser intensity [201 EI] • The characteristic energy scale then is given by the single-photon 
recoil energy 

h 2 kl . . 

Erec = ^ , (1.23) 

where M denotes the atomic mass. To give a numerical example, when working with 87 Rb 
in a lattice generated by laser radiation with wavelength A = 2-k /k^ = 842 nm [H [5] one has 
E Tec = 1.34 • 10 -11 eV. Thus, typical lattice depths of 5 to 10 recoil energies are on the order 
of 10 -10 eV — which means that one encounters many phenomena with ultracold atoms 
in optical lattices which are known from traditional solid-state physics, but scaled down in 
energy by no less than 10 orders of magnitude. 

This also tells us what "ultracold" means. Taking an ensemble of atoms with a temper- 
ature T ens such that k-oT ens is roughly equal to E rec , say, where k-Q is Boltzmann's constant, 
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the de Broglie wavelength of these atoms, given by 



h 



_2_X 
7^ 2 



(1.24) 



deBrogli 



y/2nMk B T t 



ens 



is barely longer than the lattice constant a = A/2. But in order to experience quantum 
mechanical lattice effects, the particles have to be able to "feel" the periodic structure, so 
that AdeBmgiie should cover at least a few lattice constants — which means that being this 
cold is not cold enough: We even require ^B^ens "C E rcc . 

With hardly any thermal excitation energy left the atoms occupy only the lowest Bloch 
band of their optical lattice, so that the single-particle Hamiltonian with the lattice poten- 
tial (ll.22p translates directly into the single-band tight-binding model (11 .ip when working 
in a basis of Wannier functions pertaining to that lowest band, and neglecting all couplings 
other than those between nearest neighbors, denoted as J. The accuracy of this approxima- 
tion increases with increasing lattice depth [U [T2]: For Vo/E Tec = 5 the magnitude of the 
ratio of the neglected matrix element connecting next-to-nearest neighbors to J still reaches 
about 5%, but it decreases to about 1% when V /E Tec = 10. Moreover, when expressing 
the exact band structure of a cosine lattice in terms of characteristic values of the Mathieu 
equation, and noting that the width W of the cosine energy band (jl.3p is 4J, one finds the 
approximation [21] 



The requisite still missing now is the time-periodic force corresponding to the model (II. 51) . 
This can be effectuated either by introducing a small oscillating frequency difference between 
the two lattice-generating laser beams, as detailed later, or by retro-reflecting one such beam 
off an oscillating mirror back into itself [31 IH |5] . In a frame of reference co-moving with the 
oscillating lattice, one then obtains the single-particle Hamiltonian 





(1.25) 



Hit) 




+ Vi at (x) - FtX COs(ut + 0) , 



(1.26) 



where p is the atomic center-of-mass momentum in the lattice direction, the driving force is 
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Figure 1.1: Above: One Brillouin zone of quasienergies for the optical lattice (ll.22p with 
depth V /E Tec = 5.7, driven with scaled frequency hu/E TCC = 0.5, vs. the scaled driving ampli- 
tude K . The lower left panel testifies that the first band collapse is almost perfect, whereas 
the second one, enlarged in the lower right panel, is already thwarted by multiphoton-like 
resonances. 



parametrized in accordance with Eq. (11. 9p . and we have also admitted an arbitrary phase 4>. 

In all our model calculations we consider a lattice with depth V /E rec = 5.7, implying 
that the width of the lowest Bloch band is W/E Tec = 0.22, whereas the gap between this 
lowest band and the first excited one figures as A/E rec = 2.76. Even for such a com- 
paratively shallow lattice, which is routinely being realized in current experiments [5], the 
dispersion of the lowest band already is reasonably well described by the tight-binding co- 
sine approximation (11.31) . setting J = W/4. In order to obtain dynamic localization, the 
driving frequency should then be chosen such that the quantum hu is significantly smaller 
than the gap A, so that, perturbatively speaking, interband transitions require higher-order 
multiphoton-like processes, which would be suppressed as long as the driving amplitude F\ 
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is not too strong [19] . On the other hand, it is reasonable to demand that hcu be larger than 
the band width, so that the band fits into a single quasienergy Brillouin zone. A good choice 
of the driving frequency should therefore adhere to the chain 4J = W < Hu < A; we take 
hu/E TCC = 0.5 in all numerical scenarios depicted below. For 87 Rb atoms in a lattice with 
A = 842 nm this choice fixes the frequency at u/(2tt) = 1.62 kHz. 

Figure 11.11 shows one Brillouin zone of quasienergies for these parameters vs. the scaled 
driving amplitude Kq, as defined by Eq. (I1.16p . Observe that the first quasimomentum 
Brillouin zone ranges from —hit /a = —hk^ to -\-hir/a = +hk^; the homogeneous force does 
not mix states with different wave numbers [10] . Hence, we combine quasienergies for states 
with k = (z/10)/cl in this plot, with i = 0, 1, 2, . . . , 10. In this way, the comparison of 
the ideal quasienergy band (11. lip with the one appearing in the actual optical lattice is 
greatly facilitated. Evidently the first band collapse is almost perfect, although it is slightly 
shifted from K = 2.405, the first zero of J , to K m 2.35. In contrast, the second collapse, 
expected at K = 5.520, already is significantly affected by a host of anticrossings, indicating 
multiphoton-like resonances. Thus, with V /E Tec = 5.7 and hu/E rcc = 0.5 we may expect 
almost perfect dynamic localization at the first collapse point, whereas there will be strong 
disturbances of the ideal dynamics at the second one. 

In Fig. I1.2l we depict the lowest quasienergy band for K = 0, where it coincides with the 
original energy band; K Q = 1.18, where its width is reduced by a factor of J (1.18) = 0.681; 
and at the first collapse point, K = 2.35. Ideally, a collapsed quasienergy band is completely 
flat, so that dynamic localization is associated with an infinite effective mass of the driven 
Bloch particle. Here we still observe some residual dispersion, probably resulting from both 
next-to-nearest neighbor couplings and couplings to higher bands, but the degree of band 
flattening achieved by the driving force is nonetheless impressive. 

The ultimate demonstration of dynamic localization requires, of course, the inspection of 
wave-packet dynamics. To this end, we first compute the Bloch states (xlfx^) of the lowest 
energy band of the lattice (ll.22p . and use them to design an initial wave packet 




(1.27) 
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Figure 1.2: "Lowest" quasienergy band for the optical lattice ( jl,22p with depth V / 'E rec = 5.7, 
driven with scaled frequency hu/E rec = 0.5, and scaled amplitudes K = (left), 1.18 
(middle), and 2.35 (right). Additional curves result from higher bands. 



with a Gaussian fc-space distribution 



centered around some predetermined wave number k c , with width Ak. The correspond- 
ing probability density \{x\ip{t = 0))| 2 is concentrated in the wells of the lattice potential, 
equipped with a Gaussian envelope that varies the more slowly with x the narrower its dis- 
tribution (11.281) . that is, the smaller Ak. We then take this packet (I1.27P as initial condition, 
and compute the wave function (x\ip(t)) for t > by solving the time-dependent Schrodinger 
equation numerically, fixing the phase <fi in the Hamiltonian (I1.26P at the value <fi = tt/2. 
This means that the force F(t) = F\ cos (cut + <j)) is instantaneously switched on at t = 0. 
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z/jt 

Figure 1.3: Spreading of the Bloch wave packet (I1.27P with initial fc-space width Ak/k^ = 0.1, 
and initial momentum k c /k^ = 0, in the unforced optical lattice. In this and the following 
figures, density is encoded in shades of gray. 




Figure 1.4: Evolution of the same initial wave packet as in Fig. 11.31 at the first band col- 
lapse (Kq = 2.35): Here one encounters almost perfect dynamic localization; wave-packet 
spreading is disabled because the quasienergy band is dispersionless. 
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Figure 1.5: Evolution of the same inital wave packet as in Fig. 11.31 at the supp osed second 
band collapse (K Q = 5.52): Here the multiphoton-like resonances visible in Fig. 11.11 lead to 
a marked degradation of the localization. 

Figure [L3l shows a density plot of the wave packet when it evolves in the undriven lattice, 
that is, for Kq = 0; the density is encoded in shades of gray. In this and the following figures, 
spatial extensions are measured in terms of the dimensionless coordinate z = k^x, so that a 
distance Az/tx = 1 corresponds to one lattice period; moreover, the time scale is set by the 
period T = 2n/u. With k c /k-^ = the initial packet carries no net momentum; its width is 
chosen as Ak/k-^ = 0.1. As expected, the width of the packet then grows in the course of 
time by well-to-well tunneling. 

In Fig. 11.41 we depict the density of the wave packet that evolves from the same initial 
condition when the driving amplitude is tuned to the first band collapse at K = 2.35. Here 
we observe dynamic localization at its very best: The spreading has stopped, the packet is 
"frozen" . 

It is then also of interest to monitor the evolution at the supposed second collapse, at 
Kq = 5.52; this is done in Fig. 11.51 While the "regular spreading" that has been prominent in 
Fig. ll.3l indeed seems to have stopped, small probability wavelets leak out of the initial packet 
almost immediately, spreading rapidly over the lattice. This is an effect of the multiphoton- 
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like resonances previously spotted in Fig. II. 1[ which assist parts of the wave function in 
getting to higher bands, allowing them to escape on a short time scale. 

As long as interband transitions remain negligible, the resulting single-band dynamics can 
often be regarded as "semiclassical" [22]: Namely, if an initial packet is strongly centered in 
fc-space around some arbitrary wave number k c = k c (0), this center wave number evolves in 
time according to Bloch's famous "acceleration theorem" 

hk c (t) = F(t) , (1.29) 

similar to the evolution ( II. 7p of the index of a single Houston state. The model Hamilto- 
nian (11.261) specifies F(t) = Fi cos(cot + 0), so that in this case 

k c (t) = fc c (0) + ^(sin(u;t + 0) - sin(0)) . (1.30) 

The packet's group velocity then is given by the derivative of the dispersion relation E(k) 
of the band it lives in, evaluated at this moving center wave number f 1 1 . 3 j) : 



ld£ 



(1.31) 

fcc(t) 



Taking the tight-binding relation ( 11. 3D as a good approximation for the actual lowest energy 
band of our model, this yields 

2Ja 



(t) = — — sin(fc c (t)a) . (1.32) 



^group v " ' * 



Upon time-averaging, one is therefore left with 

_ 2J c ga 



h 



sin (k c (0)a- K sm((j)) ), (1.33) 



where J e g again is the driving-dependent effective hopping matrix element ( 11.121) . and Kq is 
the scaled amplitude (11.161) . Thus, the initial phase (f) may be utilized for imparting some 
momentum to the packet. Nonetheless, for any combination of k c (0) and <fi the average group 
velocity vanishes when J e fj = 0, as corresponding to ideal dynamic localization. 
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Figure 1.6: Evolution of the fc-space distribution initially given by Eq. (jl.28p with width 
Ak/k-L = 0.1, for Kq — 1.2 and k c /k^ = 0. The white-dashed line indicates the "classical" 
solution fOOl) . 



5 








Figure 1.7: Evolution of the fc-space distribution initially given by Eq. ( 11.28P with width 
Ak/k^ = 0.1, but now for K = 0.4 and k c /k-^ = 0.8, so that the white-dashed "classical" 
solution fll.30p starts from a non-zero value. 
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2.1% 



Figure 1.8: Evolution of the Bloch wave packet (11.2 T[) with initial /c-space width Ah/k^ = 0.1 
and initial momentum /c c //cl = 0.8, driven with scaled amplitude K = 0.4, as corresponding 
to the fc-space distribution depicted in Fig. 11.71 Th e white line marks the trajectory obtained 
by integrating the oscillating group velocity (11.32D . 




Figure 1.9: Evolution of the Bloch wave packet ( 11. 27ft with initial £;-space width Ak/k-^ = 0.1 
and initial momentum k^jk^ = 0.8, now driven with scaled amplitude Kq = 2.35: Despite 
the nonzero average momentum, the average group velocity vanishe s be cause of the band 
collapse. The white line in the center is obtained as described in Fig. 11.81 
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This semiclassical behavior is illustrated by a further set of figures. In Fig. 11.61 we plot 
the evolution of the exact /c-space density that originates from the initial condition (jl.28p . 
Again we set Ak/k^ = 0.1, meaning that the distribution is sufficiently narrow to ensure the 
validity of Eq. ( jl,3ip ; moreover, k c /k^ = and K = 1.2. Since = 7r/2, the distribution 
then oscillates around k = —Fi/(hu>), or k/k^ = —K /tt, following precisely the fc-space 
trajectory predicted by Eq. (11.301) . 

A nonzero average momentum of the packet can likewise be achieved by selecting some 
suitable value of k c /k^. Figure fl~7l shows an example with k c /ki, = 0.8, while Kq = 0.4 and 
Ak/k-L = 0.1. This obviously corresponds to a wave function (x\tp(t)) which moves into the 
positive x-direction all the time; the density of this wave function is displayed in Fig. I1.81 
Here the white line indicates the classical trajectory that results from integrating the group 
velocity (11.320 ; indeed, this trajectory describes the motion of the packet's center quite well. 
When adjusting the driving amplitude to the first collapse, as in Fig. 11.91 the average motion 
stops despite the nonzero average momentum, as it should; when increasing K to still higher 
values, so that J e g becomes negative, the packet's direction of motion can even be reversed. 

While the semiclassical approach to dynamic localization may be helpful, insofar as it 
appeals to our intuition, its explanation in terms of "prohibited dephasing" resulting from a 
quasienergy band collapse is much more powerful: This view immediately reveals that not 
only does the average motion of a wave packet come to a complete standstill, but so does 
its spreading; moreover, prohibited dephasing applies to any initial condition, regardless 
whether or not its envelope varies suffiently slowly to justify the semiclassical approxima- 
tion. As an extreme example of "nonclassical" motion we consider in Fig. 11.101 the undriven 
evolution of a wave function that coincides with a single Wannier function of the optical 
lattice [12] at t — 0, and therefore certainly does not possess a slowly varying envelope then, 
giving rise to a fairly complex spreading pattern which differs substantially from the semi- 
classical one previously visualized in Fig. 11.31 Nonetheless, when driven with the amplitude 
Kq = 2.35 marking the first quasienergy band collapse, one observes another occurrence of 
dynamic localization, as witnessed by Fig. 11.111 the difference between the two evolution 
patterns depicted in Figs. 11.101 and 11.111 could hardly be more striking. 
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Figure 1.10: Evolution of the wave function that originates from an initial single Wannier 
state in the undriven lattice. This state does not possess a slowly varying envelope, and thus 
does not conform to semiclassical dynamics. 




z/% 

Figure 1.11: Evolution of the wave function that originates from an initial single Wannier 
state when driven with scaled amplitude Kq = 2.35. The semiclassical approximation is not 
applicable here, but dynamic localization works nonetheless. 
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Figure 1.12: Experimental setup for in-situ measurement of dynamic localization of a Bose- 
Einstein condensate (BEC) in a driven optical lattice: The frequencies of two laser beams 
are shifted with the help of acousto-optic modulators (AOMs) by v and by v + Au sm(ut) , 
respectively, before being directed against each other by mirrors. The resulting optical lattice 
then oscillates in the laboratory frame, giving rise to an oscillating inertial force in the frame 
of reference co-moving with the lattice. After the initial longitudinal confinement is switched 
off, the BEC expands by well-to-well tunneling; its final width (indicated by the dashed line) 
is recorded by imaging its shadow cast by a resonant flash onto a CCD chip. (Figure courtesy 
of O. Morsch.) 



In actual laboratory experiments it is advantageous to work with a phase-coherent atomic 
Bose-Einstein condensate, rather than with individual atoms: If the density of the condensate 
is sufficiently low, or if the interatomic s-wave scattering length is tuned close to zero by 
means of a Feshbach resonance [13], the condensate is practically ideal, so that one effectively 
can perform a measurement on an ensemble of identically prepared noninteracting atoms in 
a single shot. Figure [LT2l shows a possible experimental setup [3j H]: The optical lattice 
is formed by two laser beams of wavelength A, which are directed against each other with 
the help of mirrors. Each beam passes through an acousto-optic modulator which shifts its 
frequency by v and by u+ Ai>(t), respectively. Because of the frequency difference Av(t) thus 
introduced between the counterpropagating beams, the condensate experiences the potential 



V ( 



A /■* 

x + - / dr Avir) 
I Jo 



1.34) 



in the laboratory frame, which means that the lattice position shifts in time according to 
the prescribed protocol Au(t). In a frame of reference co-moving with the lattice, this shift 
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translates into the inertial force 

f( t) = M^^W. (1 .3B) 
w 2 dt K ' 

Therefore, choosing Au(t) = Ai/ max sin(o;t + 0) leads to the desired Hamiltonian (11.261) in 
the co-moving frame, with the driving amplitude 

Fi = Mw ^ Az/ max . (1.36) 

Now a Bose-Einstein condensate initially trapped in the center of the oscillating lattice is 
allowed to expand freely in the lattice direction by well-to-well tunneling after switching off 
the longitudinal confinement, while maintaining a weak transversal confinement in order to 
keep the condensate in the lattice. After a variable expansion time, the in situ width of the 
condensate is determined by a resonant flash, the shadow cast by which is imaged onto a CCD 
chip |3j. The measured expansion rate then is to a good approximation proportional to | J e ff| ; 
that is, to the absolute value of the effective hopping matrix element (11.121) : in principle, 
even the sign of J c g can be deduced from additional time-of-flight measurements |3j. In 
Fig. 11.131 we display data for the ratio J e s/J acquired in this manner by the Pisa group 
with a condensate of 87 Rb atoms in a lattice of depth Vq/ E rec = 6.0, driven with frequency 
bjj (2ir) = 4.0 kHz, after expansion times of 150 milliseconds. Evidently these data match the 
expected Bessel function Jo(-Ko) quite well even up to the second zero. Note that here one 
has hu/E iec = 1.24, so that the frequency employed in these measurements is significantly 
higher than in our model calculations. This means that the inequality 4J < hu is satisfied 
in a stronger manner, while hu still remains reasonably small compared to the band gap. As 
a consequence, even the second band collapse can be quite well developed. In any case, this 
figure strikingly demonstrates that the concept of dynamic localization by now has crossed, in 
the context of mesoscopic matter waves, the threshold from an idealized theoretical concept 
to a well-controllable laboratory reality. 
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Figur e 1.13 : Experimental results for the ratio J e s/J of the effective hopp ing matrix ele- 
ment (I1.12p to the bare one as a function of the scaled driving amplitude (I1.16p . obtained with 
a Bose-Einstein condensate of 87 Rb atoms in an optical lattice with a depth of Vq = 6 E Tec 
(A = 842 nm), driven with frequency u/(2n) = 4.0 kHz. The dashed line corresponds to the 
expected Bessel function J (K ). (Figure courtesy of O. Morsch.) 
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1.3 What is it good for? 



Up to this point we have considered no more than a possible realization of dynamic local- 
ization which comes fairly close to the theoretical ideal pp. Apart from its observation with 
dilute Bose-Einstein condensates in time-periodically shifted optical lattices [31 H] , this type 
of quantum wave propagation has meanwhile also been made visible by means of an optical 
analog based on sinusoidally-curved lithium-niobate waveguide arrays [23] . This is certainly 
interesting, but it is not what one would call "deep"; the "prohibited dephasing"-view clearly 
reveals that the only physics entering here is summarized by stating that an initial state is 
"frozen" in time if the phase factors of all of its spectral components evolve at the same 
speed. Yet, the accompanying band collapse furnishes a strong hint that there may be more 
in stock. Namely, when the ideal dynamics is somehow perturbed it is the bandwidth which 
sets the scale with respect to which the strength of such a perturbation has to be gauged. A 
prominent example is provided by the repulsive interaction between ultracold atoms in an 
optical lattice; the strength of this interaction is expressed in terms of a parameter U which 
quantifies the repulsion energy of one pair of atoms occupying the same lattice site (21] • Ac- 
cordingly, the characteristic dimensionless parameter then is the ratio U/ J; here J = W/ 4 
is taken instead of the bandwidth W. Indeed, it is this ratio U/J which decides which quan- 
tum phase a gas of ultracold, repulsively interacting atoms in an optical lattice adopts: For 
U/J< 1 the system is superfluid, but becomes a Mott insulator when this ratio exceeds a 
critical value [21]. Hence, when recalling that J is replaced by the effective hopping matrix 
element f II . 1 2 j) when the system is driven with appropriate parameters, it is only natural to 
predict that this superfluid-to-Mott insulator transition can be induced in a time-periodically 
shifted optical lattice by varying the driving force [2J1 [25] : Assuming that one starts in the 
superfluid phase, J e g can then virtually be made arbitrarily small by adjusting the scaled 
amplitude Kq to a zero of Jo, resulting in a value of U/J e g so large that the system is forced 
to enter the Mott regime. The experimental confirmation of this scenario, achieved by the 
Pisa group [5], probably constitutes the first known example of coherent control exerted by 
means of time-periodic forcing on a quantum phase transition. 

There are other types of perturbations, associated with deviations from perfect transla- 



1.3. WHAT IS IT GOOD FOR? 



23 



tional symmetry, which affect even noninteracting ultracold atoms in optical lattices. Most 
notably, the system governed by the tight-binding Hamiltonian 



differing from its antecedent (11. ip through additional on-site energies which oscillate along 
the lattice with amplitude V, shows a quite peculiar behavior when the number g is irrational, 
so that this system becomes quasiperiodic [261 [28] : As long as \V/J\ < 2, so that the 
on-site perturbations are relatively weak, all of its energy eigenstates still remain extended 
over the entire lattice in a Bloch-like manner, whereas they are all exponentially localized, 
with one common localization length, when \V/J\ > 2. Thus, there is a metal-insulator-like, 
incommensurability-induced transition at \V/J\ = 2, originally studied by Harper [26] in 
the context of conduction electrons in a magnetic field, and later by Aubry and Andre |27J; 
this transition can be realized approximately with ultracold atoms in a bichromatic optical 
lattice described by the potential 



The guiding idea here is to employ a primary lattice with depth Vq for setting up the hosting 
tight-binding system (11. ip . as before, and then to invoke a secondary lattice with much 
smaller depth 1V\ for achieving the required modulation of the local energies at the sites of 
the host [TOJQT]. When the primary lattice is comparatively shallow, possessing a depth of 
only a few recoil energies, the transition occurs stepwise upon increasing V\ [12], featuring 
pronounced mobility edges resulting mainly from the next-to-nearest neighbor couplings 
between the host's sites which are present in the full bichromatic potential (11.381) . but do 
not occur in the Aubry- Andre model fll.37]) . When V /E rec ^> 1, so that the primary lattice 
is so deep that these additional couplings may be safely neglected, the transition occurring 
in the actual bichromatic lattice (11.381) is fairly sharp. The parameter J then again is given 
approximately by Eq. ( 1 1.251) : moreover, one has 



Haa = — 



jJ2(\£ + l)(£\ + \£)(£+l\) + VY,cos(2ng£ + 5)\e)(e\ , 



(1.37) 




Vbict^) — tt cos (2k-Lx) + Vi cos(2gk^x + 5) . 



(1.38) 





(1.39) 
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with reasonably chosen g on the order of unity. Therefore, the equation \ V/J\ = 2 marking 
the metal-insulator-like transition in the ideal Aubry-Andre model now translates into the 
estimate [12] 

»l „ » (K) 3, \J- 2 .rf + ? \ (1 . 40) 

E rec VVT \E Iec J y V ^rcc yVo/£ rec / 

for the critical strength of the secondary optical lattice, given a sufficient depth of the 
primary one. Indeed, this transiton has been observed with a Bose-Einstein condensate 
consisting of 39 K atoms, using a magnetically tunable Fesbach resonance for rendering these 
atoms practically noninteracting [T3] . 

When ultracold atoms in such a bichromatic lattice dl .38[) are subjected to time-periodic 
forcing, one obtains an additional knob which can be turned to induce the transition: Be- 
cause J is replaced by the effective hopping strength (II . 12[) when the system is suitably 
driven, one can cross the critical border \V/J^[\ = 2 by varying the parameters of the driv- 
ing force; the critical parameters then are linked approximately by the relation 

Jo ^0 S — ^ — 1 — ) eXP +2 V P / ; • L41 

8 E vcc \E tec J y V E mc ^V /E iec J 

Hence, it is feasible to coherently control the metal-insulator-like transition exhibited by 
noninteracting ultracold atoms in properly designed bichromatic optical potentials through 
time-periodic forcing (TOj [TT]. In order to substantiate this prediction, we now display the 
results of further numerical wave-packet calculations. In all of these we employ a primary 
lattice with depth Vq/ 'E rec = 5.7, as in our preceding studies, and fix the incommensurability 
parameter at the golden mean g = (v^5 — l)/2 up to numerical accuracy. With this choice, 
the above estimate f ll.40p yields Vi/E rec w 0.165 for the critical strength of the secondary 
lattice. The driving frequency is given by hu/E rcc = 0.5 throughout. 

Figure 11.141 visualizes the evolution of a wave function that originates from the same 
Gaussian initial state as already employed in Fig. 11.31 Here the driving force is still absent, 
and the depth of the secondary lattice is V\/E rec = 0.10, placing the system in its metallic 
phase; accordingly, the wave function readily explores the entire lattice. In contrast, when 
V\/E xec = 0.25 and the drive is still switched off, the wave function remains localized as shown 
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Figure 1.14: Evolutio n of t he same initial wave packet as in Fig. 11.31 in an undriven trichro- 
matic optical lattice (|1.38p . Here the strength of the secondary potential is Vi/E rec = 0.10, 
so that the system is in its mobile "metallic" phase, allowing the wave function to spread. 




Figure 1.15: Evolutio n of t he same initial wave packet as in Fig. 11.31 in an undriven bichro- 
matic optical lattice ( 11.38R . Here the strength of the secondary potential is V\/E iec = 0.25, 
so that the system is in its "insulating" phase, keeping the wave function localized. 
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Figure 1.16: Evolution of the same initial wave packet as in Fig. 11.31 in a Dichromatic optical 
lattice (jl.38p driven with s caled amplitude K = 1.7. The strength of the secondary lattice 
is V\/E Tec = 0.10, as in Fig. 11.141 so that the system would be in its "metallic" phase if there 
were no forcing. 

in Fig. I1.15| this indicates that we are encountering the insulating phase now. But the wave 
function also remains localized when the secondary lattice is tuned back to Vi/E Iec = 0.10 
and the driving force acts with scaled amplitude K = 1.7, as depicted in Fig. 11.161 The 
relation (11.411) predicts the transition from the metallic to the insulating phase to have 
occurred already at about K « 1.3. It should be noted that there is a pronounced difference 
from the ideal dynamic localization reviewed in the preceding section: There the wave packet 
remains localized only when Kq is exactly equal to a zero of Jo- In contrast, here one switches 
from the metallic into the insulating phase already when |Jo(i£o)| becomes sufficiently small. 

Finally, we show a corresponding sequence of results for wave functions which evolve from 
an initial Wannier state of the primary lattice. In Fig. 11.171 we again consider an undriven 
bichromatic lattice with V\/E XfX = 0.10, so that the mobile metallic phase enables uninhibited 
spreading; in Fig. 11.181 where Vi/E vec = 0.25, the system's insulating character then keeps 
the wave function strongly localized. But that same high degree of localization may also be 
obtained when again resetting the strength of the secondary lattice to V\/E TCC = 0.10, and 
switching on the driving force with scaled amplitude Kq = 1.7, as done in Fig. 11.191 
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Figure 1.17: Evolution of the wave function originating fro m a single Wannier state of the 
primary lattice in an undriven bichromatic optical lattice (ll.38p . Here the strength of the 
secondary potential is V\/E rec = 0.10, so that the system is in its "metallic" phase. 
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Figure 1.18: Evolution of the wave function originating fro m a single Wannier state of the 
primary lattice in an undriven bichromatic optical lattice (II. 38ft . Here the strength of the 
secondary potential is V\/E iec = 0.25, so that the system is in its "insulating" phase. 
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Figure 1.19: Evolution of the wave function o rigina ting from a single Wannier state of the 
primary lattice in a bichromatic optical lattice ( jl.38j) driven with sc aled a mplitude K = 1.7. 
The strength of the secondary lattice is Vx/E rec = 0.10, as in Fig. 11.171 so that the system 
would be in its "metallic" phase if there were no forcing. 

These figures vividly illustrate the main message: In the presence of time-periodic forcing 
it is the width of the underlying quasienergy band which determines the effective strength 
of deviations from perfect spatial periodicity. In an ideal lattice without such deviations 
one encounters "only" dynamic localization, but in lattices with isolated, quasiperiodic, or 
random perturbations the strengths of these can be adjusted at will by suitably selecting 
the parameters of the drive. With regard to experimental tests, the enormous flexibility 
offered by ultracold atoms in optical potentials makes such systems far superior to electrons 
in ac-driven crystal lattices. 

When the concept of controlling the incommensurability-induced metal-insulator transi- 
tion exhibited by the Aubry-Andre model (11.371) by means of time-periodic forcing was con- 
ceived [TQl [11] the experimental investigation of ultracold atoms in optical lattices was still 
in its infancies. But now that this transition has been unambiguously observed with a non- 
interacting Bose-Einstein condensate [13], the demonstration of its coherent control has come 
into immediate reach. Besides the already established coherent control of the interaction- 
induced superfluid-to-Mott insulator transition [5], this demonstration would constitute a 
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further milestone achievement in the on-going effort to explore the newly emerging prospects 
provided by dressed matter waves. 
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